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Abstract 

We find that localized quantum iV-body soliton states exist for a derivative non- 
linear Schrodinger (DNLS) model within an extended range of coupling constant 
given by < \£ q \ < ^ tan (7^1) ■ We a ^ so observe that soliton states with both positive 
and negative momentum can appear for a fixed value of £ 9 . Thus the chirality prop- 
erty of classical DNLS solitons is not preserved at the quantum level. Furthermore, 
it is found that the solitons with positive (negative) chirality have positive (negative) 
binding energy. 
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1 Introduction 



Quantum soliton states associated with integrable field models in 1+1 dimension have 
attracted a lot of interest in recent years [1-10]. These quantum soliton states are 
usually constructed by using the approach of either the coordinate Bethe ansatz or 
the algebraic Bethe ansatz. By applying the coordinate Bethe ansatz to an integrable 
nonrelativistic Hamiltonian, one can directly obtain its exact eigenf unctions in the 
coordinate representation. If such an eigenfunction decays sufficiently fast when a rel- 
ative coordinate measuring interparticle distance tends towards infinity, we define this 
localized square-integrable eigenfunction as a quantum soliton state. It is also possi- 
ble to construct quantum soliton states through the approach of the algebraic Bethe 
ansatz, by choosing appropriate distributions of the spectral parameter on the complex 
plane [2-4]. The stability of quantum soliton states, in the presence of small external 
perturbations, can be measured by calculating their binding energy. It is usually found 
that localized quantum soliton states of various integrable models, including the well 
known nonlinear Schrodinger model and sine-Gordon model, have positive binding en- 



The purpose of the present article is to study the quantum soliton states of an 
integrable derivative nonlinear Schrodinger (DNLS) model [6-10]. Classical and quan- 
tum versions of DNLS model have found applications in different areas of physics like 
circularly polarized nonlinear Alfven waves in a plasma [11,12], quantum properties of 
the optical solitons in fibers [13] and some chiral Luttinger liquids which are obtained 
by dimensional reduction of a Chern-Simons model defined in two dimensions [14-15]. 
It is well known that the classical DNLS model can have solitons with momenta in only 
one direction [16-18]. Here we want to investigate whether this chirality property of 
the classical solitons is preserved at the quantum level. The Hamiltonian of quantum 
DNLS model in its second quantized form is given by [6,7] 



where we have set the particle mass m = 1/2, £ 9 is the coupling constant, and 
ip{x), ijj'(x) are field operators obeying the following equal time commutation rela- 
tions: [ip(x),ip(y)} = [il)'(x),i/)'(y)] = 0, [i[>(x),i/)'(y)] = h8(x — y). Quantum soliton 
states of this DNLS model have been constructed through the approach of the alge- 
braic Bethe ansatz [8-10] as well as the coordinate Bethe ansatz [6,7]. By applying the 
coordinate Bethe ansatz, it is found that quantum iV-body soliton (henceforth called 
iV-soliton) states exist for this DNLS model provided that £ q is restricted within the 



ergy [1-5]. 




dx \ ipl(x)ip x (x) +i^ q {^ 2 (x)^(x)^ x (x) - ipl(x)^(x)ip 2 (x)} 1 , (1.1) 
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range: < \£ q \ < |tan^j. Moreover it is observed that, similar to the classical 
case, such iV-soliton states can have only positive momentum for the case £ 9 > [6]. 
But, it was mentioned later that, soliton states can exist even for negative momentum 
for small values of N [7]. However, no explicit range of £ 9 is known for which such 
7V-solitons with negative momenta may appear in the quantum DNLS model. 

In this article we reinvestigate the range of £ g for which localized quantum N- 
soliton states exist in the case of DNLS model. In Sec. 2, we apply the coordinate 
Bethe ansatz for this purpose and show that quantum iV-soliton states exist within an 
extended range given by < < ^tan (]y~r)- We also find that the momentum 
of the iV-soliton state has the same sign as the coupling constant when it lies within 
the range < \£ q \ < |tan(-^) and has the opposite sign to the coupling constant 
when |tan(j0 < < |tan(-^-j-). As a result, soliton states with both positive 
and negative momentum can appear for a fixed value of £ q . In Sec. 3 we find out the 
binding energy of these soliton states. Sec. 4 is the concluding section. 



2 Quantum N-soliton states for DNLS model 

For applying the coordinate Bethe ansatz, we separate the full bosonic Fock space asso- 
ciated with Hamiltonian (1.1) into disjoint N-particle subspaces \Fn). We are interested 
in solving the eigenvalue equation H\Fn) = E\Fn). The coordinate representation of 
this equation yields 

H N t n (x 1 , x 2 ,---,x N ) = E r N (x 1 , x 2 , • ■ ■ , x N ) , (2.1) 

where the N-particle symmetric wave function tn(xi, x 2 , • • • , xn) is defined as 

tn(xi,x 2 , ■ ■ • ,x N ) = l/v / n!(0|^(xi) • • • ip(x N )\F N ) , (2.2) 

and if at, the projection of the second-quantized Hamiltonian H (1.1) on the N particle 
sector, is given by 

Hn = ~h 2 £ ^ + 2ih% £ S( Xl - x m ) (— + — ) . (2.3) 

3=1 UX j Km UXl UXm 

It is evident that this H N commutes with the total momentum operator in the iV- 
particle sector, which is defined as 

N Q 
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Next, we divide the coordinate space R N = {xi,x 2 , ■ ■ -xn} into various iV-dimensional 
sectors defined through inequalities like < 2^(2) < • • • < £ w (iv)> where cj(1),u;(2), 
■ ■ ■ , u(N) represents a permutation of 1, 2, • • • , N. Within each such sector the inter- 
action potential of the Hamiltonian (2.3) is zero and the resulting eigenfunction is just 
a superposition of the free particle wave functions. Coefficients associated with these 
free particle wave functions can be obtained through the interaction part of the Hamil- 
tonian (2.3), which is nontrivial at the boundary of two adjacent sectors. It is known 
that all such necessary coefficients for the Bethe ansatz solution of a iV-particle system 
can be obtained by solving the corresponding two particle problem [19]. Hence let us 
first try to construct the eigenf unctions of the Hamiltonian (2.3) for the two particle 
case, without imposing any symmetry property on t 2 (xi,x 2 ) under the exchange of 
particle coordinates. For x\ < x 2 region, we may take such eigenfunction as 

r 2 (x 1 ,x 2 ) = exp {i(kiXi + k 2 x 2 )} (2.5) 

where k\ and k 2 are two distinct wave numbers. Using eqn. (2.1) for N = 2, one can 
show that this two-particle wave function will take the following form in the region 

xi > x 2 : 

t 2 (xi,x 2 ) = A(k u k 2 ) exp {1(^X1 + k 2 x 2 )} + B(k 1 , k 2 ) exp{i(k 2 x l + hx 2 )} , (2.6) 
where the 'matching coefficients' A(k±, k 2 ) and B(k±, k 2 ) are given by 

A(k 1 ,k 2 )= ih - h ]^ lh lf 1 + h) , B(k l ,k 2 ) = l-A(k l ,k 2 ). (2.7) 
(fci - k 2 ) 

By using these matching coefficients, one can construct completely symmetric N- 
particle eigenfunctions for the Hamiltonian (2.3). In the region x\ < x 2 < ■ ■ ■ < xn, 
such eigenfunctions are given by [6,19] 

t n (x 1 ,x 2 ,---,xn) =Y1 f II —J^T^ 1 P^)M2l-MN)(xu%2r ■ ■ ,%n) , (2.8) 

w \l<m A K K m, K l) J 

where 

Po;(l),o;(2),-,a;(Ar)(^l,^2, • • - ,X N ) = exp H h ^(N^n)} ■ (2.9) 

In the expression (2.8), k n s are all distinct wave numbers and J2u implies summing over 
all permutations of the numbers (1,2,....N). The eigenvalues of the momentum (2.4) and 
Hamiltonian (2.3) operators, corresponding to the eigenfunctions tn(xi, x 2 , • ■ ■ , %), 
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are given by 

N 

P N t n (x 1 , x 2 , ■ ■ ■ , x N ) = h(Yl k i) t n{xi,x 2 , ■ ■ ■ , x N ) , 

3=1 
N 

H N r N (x 1 ,x 2 , ■ ■ -,x N ) = h 2 (J2 k ]) t n(xi,x 2 , ■ ■ -,x N ) . (2.10a, b) 

3=1 

The wave function (2.8) will represent a localized soliton state, if it decays when any 
of the relative coordinates measuring an interparticle distance tends towards infinity. 
To obtain the condition for the existence of such localized soliton state, let us consider 
the following wave function in the region x 1 < x 2 < ■ ■ ■ < Xjy : 

N 

Pi,2,-,n{x!,x 2 , - ■ ■ ,x N ) = exp(iJ2kjXj). (2.11) 

3=1 

As before, the momentum eigenvalue corresponding to this wave function is given by 
hJ2jLikj- Since this momentum eigenvalue must be a real quantity, one gets the 
condition 

N 

5> = 0, (2.12) 

i=i 

where qj denotes the imaginary part of kj. The probability density corresponding to 
the wave function Pi,2,--,at( xi, x 2 , • • • ,xn) (2.11) can be expressed as 

N-l r 

|pi,2,...,Ar(^i,^2,---,2;iv)| 2 = exp {2 (E^)^}' ( 2 - 13 ) 

r=l j=l 

where the y r 's are the N — 1 number of relative coordinates: y r = x r+ i — x r . It is 
evident that this probability density (2.13) decays exponentially in the limit y r — > 00 
provided that all the following conditions are satisfied: 

JV-l 

gi <0, qi + q 2 <0, ,E^<°- ( 2 - 14 ) 

j'=i 

Note that the wave function (2.11) is obtained by taking oj as the identity permutation 
in (2.9). However, the full wave function (2.8) also contains terms like (2.9) with u 
representing all possible nontrivial permutations. The conditions which ensure the 
decay of such a term with nontrivial permutation, in general, contradict the conditions 
(2.14). Consequently, to have an overall decaying wave function (2.8), the coefficients 
of all terms Pu>(i),u>(2),-,u>(N)(xi, x 2 , • • • , x^) with nontrivial permutations are required to 
vanish. This requirement leads to a set of relations like 

A(k u k 2 )=0, A(k 2 ,k 3 )=0, , A(k N -!, k N ) = 0. (2.15) 



5 



Thus the simultaneous validity of the conditions (2.12), (2.14) and (2.15) ensures that 
the full wave function tn(xi,%2, • • • , %n) (2.8) would represent a localized soliton state. 
Using the conditions (2.12) and (2.15), one can derive the expression for the complex 

kfi S clS 

fc|i = xe -i(W-2n)^ (216) 

where x is an J e ^ undetermined real parameter, and is another real parameter which 
is related to the coupling constant as 

= tan" 1 (^). (2.17) 

To obtain a unique value of from the above equation, we restrict it within the 
fundamental region: — | < < |. In this context, it may be mentioned that a relation 
equivalent to (2.16) can also be obtained through the method of the algebraic Bethe 
ansatz, when quantum soliton states of DNLS model are considered [8-10]. Next, let 
us verify whether these k n s (2.16) satisfy the conditions (2.14). Summing over the 
imaginary parts of these k n s, we can express the conditions (2.14) in the form 

X ^sin{0(iV-/)}>O, (2.18) 
sin <p 

where / G [1, 2, 3, • • • , N — 1]. For the simplest case N = 2, the above condition is 
satisfied when lies within the range < < f (— | < < 0) for the choice x > 
(x < 0). Thus any value of within its fundamental region can generate a 2-soliton 
state. Next, we consider the case where N > 3. It may be noted that, sin0, sin0Z, and 
sin{0(A r — I)} are all positive (negative) when takes value within the range < < 
-^-j- (— jv^y < < 0). Hence the condition (2.18) will be satisfied for the choice x > 
(x < 0), when lies within the range < < jf-^ ( — j^zi < < 0). Due to eqn. 
(2.17), the corresponding coupling constant £ q lies within the range < £ q < -| tan -^-j- 
(— ^tan^j- < £ q < ). As a result, we get localized iV-soliton states within the 
above mentioned range of or £ ? . Earlier it has been shown that iV-soliton state can 
be obtained for the case £ g > when lies within the range < < [6] . The present 
analysis yields a wider range of for which a iV-soliton state can be constructed. 

Next, by using eqns. (2.16) and (2.10a), we obtain the momentum eigenvalue (P) 
corresponding to these soliton states as 

„ sin 6N 

P = h X — L -. 2.19 
sm0 

Thus, for any given value of 0, the momentum of such iV-soliton state is proportional 
to the value of parameter x- Due to eqn. (2.19), it is evident that P and x have the 
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same sign when lies in the range < |0| < and the opposite sign when lies in 
the range < |0| < j^zj- However, we have seen earlier that x an d £, q must have the 
same sign within the allowed range of a iV-soliton state. Consequently, the momentum 
of the iV-soliton state has the same sign as £ q when < \<f>\ < and has the opposite 
sign to £ q when < |0| < jf-^. We may call solitons with momenta having the same 
(opposite) sign as £ q as solitons with positive (negative) chirality. It may be noted that 
the range of given by -^ < |0| < jf-^ does not lie within its fundamental region for 
the case N = 2. As a result, iV-soliton states with negative chirality can appear only 
for N > 3. 

The appearance of both positive and negative momentum solitons in quantum 
DNLS model can be described most conveniently by considering a fixed positive value 
of or £ q . Let us first take a fixed value of so that ^ is an integer. In this case 
we obtain solitons with positive momentum for the choice N — 2, 3, • • • , ^ — 1 and a 
single soliton with zero momentum for N max = j. Next we consider the case where ^ 
is non-integer and denote its integral part by ^ . In this case we obtain solitons with 
positive momentum for the choice N — 2, 3, • • • , ^ and a single soliton with negative 
momentum for N max = ^ + 1. Thus we find that the chirality property of classical 
DNLS solitons is broken at the quantum level. It may be further noted that the ratio 
of solitons with positive and negative momentum is of the order of ^. Moreover, due 
to eqn. (2.17), we get — > at the classical h — >■ limit (for any fixed value of £ q ). As 
a result, the ratio of solitons with positive and negative momentum becomes infinitely 
large at the classical limit. 



3 Binding energy of a iV-soliton state 

Here we want to calculate the binding energy for the iV-soliton states of quantum DNLS 
model within the allowed range of given by < |0| < Substituting the values 

of complex k n (2.16) in (2.10b), we obtain the energy eigenvalue (E) corresponding to 
these iV-soliton states as 

sin(20) v ; 

To calculate binding energy, we assume that the momentum P (2.19) of the iV-soliton 
state is equally distributed among iV number of single-particle scattering states. The 
real wave number associated with each of these single-particle states is denoted by ko. 
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With the help of eqns. (2.10a) and (2.19), we obtain 

X sin(0AQ 
N sm <p 

Using eqn. (2.10b), one can easily calculate the total energy for N number of such 
single-particle scattering states as 

w ^9^,9 ft 2 X 2 sin 2 6N , . 

Subtracting E (3.1) from E' (3.3), we obtain the binding energy of quantum iV-soliton 

state as 

E B {4>, N)=E>-E= h2x2sin / N \P^r ~ (3.4) 

sm (ft 1 N sm (ft cos (ft > 

It may be noted that the above expression of binding energy remains invariant under the 

transformation (ft — > —(ft. Thus it is sufficient to investigate the nature of this binding 

energy only for positive values of (ft within the range < (ft < j^zi- Substituting N = 2 

to eqn. (3.4), we obtain Eb{4>-,2) = 2h 2 \ 2 sm 1 (ft. Thus we get Eb(4> : 2) > for any 

nonzero value of (ft. For N=3, eqn. (3.4) takes the form .Eb(0, 3) = 8n 3 x sin 2 0(3 — 

4 sin 2 (ft). Here we get E B ((ft, 3) > when < (ft < f and E B ((ft, 3) < when | < (ft < f . 

So from the above examples we can guess that Eb(4>,N) > when < (ft < and 

E B < when < (ft < j^zr[. By using the method of induction, we will try to prove 

that this condition is really valid for arbitrary values of N. To this end, we write 

E B ((ft,N) (3.4) in the form 

E B ^N)= h ^^Lm,N), (3.5) 
sm (ft 

where 

t( . An sin0iV cos0iV 

fW,N) = — . 3.6 

W sm (ft cos (ft 

To use the method of induction, let us first assume that /(0, N) > 0, i.e. 

sin d> N NcosSN /n n . 

sm (ft cos (ft 

when (ft lies within the range < (ft < j^. Next we consider an identity given by 

sin0(A^+l) /sin0A^\ 

— — ■ = — cos (ft + cos (ftN. 

sm0 y sm0 J 

Using the inequality (3.7) and the above identity, we obtain 

Sin0(iV + 1) >(iV + l)cos0iV, (3.8) 
sm0 
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for < (ft < fj. Next we consider another identity given by 

cosd)(N + l) sin</>sin(</>iV) 
cos <pN = — — ■ H j— — . 

COS (ft COS (ft 

Since sin0sin(0A r )/cos0 > for < (ft < the above identity yields 

,»r (iV+ l)cos0(iV + 1) . . 

(N + l)cos6N>± '- -p S 3.9 

COS0 

for < (ft < j;. Now, combining the inequalities (3.8) and (3.9), we get 

sm0(iV + l) > (iV+l)cos0(iV + l) 
sin (ft cos (ft 

for < (ft < jj. Thus by assuming that f((ft, N) > within the range < (ft < -j^, we 
are able to show that /(0, iV + 1) > within the same range of 0. By putting iV = 2 in 
eqn. (3.6), it is easy to check that f((ft, 2) > for < (ft < f . Thus by using the method 
of induction, it is established that f((ft, N) > within the range < (ft < f° r an 
iV > 3. Now, it may be noted that > for < (ft < §. Consequently, from eqn. 

(3.5) we find that the binding energy for a iV-soliton state of quantum DNLS model 
is positive when (ft lies in the range < (ft < f^. However, it may also be noted that 
s ™n!p < within the range < (ft < As a result, from eqn. (3.5) we interestingly 
find that the binding energy for a iV-soliton state is negative when (ft lies in the range 
fj < (ft < jjzj- Thus we find examples of quantum solitons in an integrable model 
where the binding energy can be negative. 



4 Concluding Remarks 

By applying the coordinate Bethe ansatz, we have investigated the range of the coupling 
constant £ q for which localized quantum iV-soliton states exist for DNLS model. It 
turns out that such quantum iV-soliton states exist within an extended range of £ q 
given by < \£ q \ < | tan ■ We also observe that the iV-soliton state has positive 

chirality when £ q lies within the range < \£ q \ < |tan^^ and negative chirality 
when |tan(j0 < \£ q \ < ^tan(-^-j-). Thus we find that the chirality property of 
classical DNLS solitons is not preserved at the quantum level. We also calculate the 
binding energy for the quantum soliton states of this model. Interestingly, it is found 
that the binding energy for a iV-soliton state is positive when < \£ q \ < | tan (^j^J 
and negative when ^tan^^ < < -r tan ■ Thus we can conclude that the 

solitons with positive chirality form bound states with positive binding energy. On the 
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other hand, solitons with negative chirality have negative binding energy and hence 
they form anti-bound states. Consequently, solitons with negative chirality would be 
unstable in the presence of external perturbations. As a future study, it might be 
interesting to calculate the decay rate of such quantum solitons by introducing small 
perturbations to the DNLS Hamiltonian. 

Recently, we have discovered additional ranges of values of for which there are 
bound and anti-bound iV-soliton states. The complete description of these ranges can 
be found using concepts from number theory such as Farey sequences and continued 
fractions [20]. 
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